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Abstract 

We introduce a novel technique to find the asymptotic time behaviour of 
deterministic systems exhibiting anomalous diffusion. The procedure is tested 
for various classes of simple but physically relevant 1-D maps and possible 
relevance of our findings for more complicated problems is briefly discussed. 



In the last few years the phenomenon of deterministic diffusion has been widely 
investigated: as a matter of fact a major theoretical challenge in the study of 
dynamical systems is to understand thoroughly the generation of typical stochas- 
tic properties in purely deterministic systems. Moreover even the simplest maps 
for which deterministic diffusion has been observed are supposed to model the 
behaviour of relevant physical systems (like Josephson junctions in presence of a 
microwave field, see e.g. U). A recent approach, introduced independently in JJ 
and [Q], leads to an expression for the diffusion coefficient in terms of the periodic 
orbits of the system, in the form of a cycle expansion^. Cycle expansions have 
been applied in a number of different contexts (see for instance |5j, |[]): they work 
remarkably well for low dimensional hyperbolic systems, provided their topology 
(symbolic dynamics) is under control. The problem of controlling the topology of 
the system is well illustrated, in the context of diffusion properties via cycle expan- 
sions, when one deals with the infinite Lorentz gas with bounded horizon (see 0). 
If one relaxes the hypothesis of pure hyperbolicity (absence of marginal stability), 
much more care has to be taken, and generally the effectiveness of the expansions 
is regained only if one is able to sum infinite contributions shadowing the marginal 
fixed point 0. We emphasize that this problem is of paramount relevance, as it 
naturally arises when dealing with generic hamiltonian systems, in which elliptic 
islands of stability and hyperbolic homoclinic webs coexist. 

Here we address explicitly the problem of marginal stability, which is tightly 
connected to the appearance of anomalous diffusion (and also represents a crucial 
feature in establishing a good theory for diffusion in generic two-dimensional area 
preserving mapsH). We extend the theory developped in 0, |3| to study deviations 
from normal diffusion, expressing the asymptotic time behaviour in terms of prop- 
erties of the periodic orbits of the system. The method is then applied to a class of 
1-D maps for which anomalous diffusion has been previously observed. 

We begin by recalling the essential features of the periodic orbit theory of normal 
diffusion ^, by considering the simplest context in which it may be applied. We 
will consider lifts of one-dimensional circle maps 

%i = /W<6N f(x + n) = n + f(x) f(-x) = -f{x) (1) 

together with corresponding torus maps x t+ \ = f(x t )\ mo di = f( x t)- Normal diffu- 
sion means that asymptotically cr 2 (t) =< (xt — xq) 2 >~ 2Dt, where the average is 
over initial conditions (by symmetry x may be taken in the unit interval). The set of 
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all periodic orbits of / will be denoted by {p}: each orbit will be characterized by its 
period n p stability A p (product of the derivatives along the cycle) and integer wind- 
ing number a p (such that for each cycle point x^ we have f np (x^ p )) = x^) + a p ). 
We focus our attention on the generating function < e^ Xt ~ x< ^ >= Q t (/3); it has been 
shown & 01 that its asymptotic behaviour is dominated by the leading eigenvalue 
of an appropriate transfer operator: Qt(/3) ~ z(0)~ t , where z(/3) is the smallest 
solution of 

c„-« = n (i - = o 

Expansion of the generating function for small j3 thus yields 

Anomalous diffusion is associated to a vanishing or diverging D. In order to 
generalize the theory we take into account that, up to time t, f2 t (/3) gets contribu- 
tions from periodic orbits with n p < £@' El (as it is connected to the trace of the 
t-th power of the transfer operator): we thus may write in general 

0*09) ~ [* w (/?)]"* (2) 

where z^(/3) is such that Co\t]( z (P) ^ P) — 0, Cow denoting the cycle expansion of 
CcT 1 truncated to finite order t. From ([|) we thus get the asymptotic behaviour 



a 2 (t) 1 



dyopH^P) 



d/dzQ^fi) 



(3) 



where we took into account probability conservation^ (z^(0) — > 1 as t — >• oo) 
as well as symmetry of the map. This procedure may be alternatively interpreted 
as follows: for a sequence of large increasing times tj. we approximate the system 
by a sequence of hyperbolic (regularly diffusing) systems with polynomial zeta 
functions Co^O, P)s k = (o[t k ]( z i P) : °" 2 (0 ~ 2D(t k ) ■ t for t ~ t k (for both the 
original system and Sj-), thus corrections to normal diffusion are connected with 
the asymptotic behaviour of D{t). 

To test the theory we consider a map of the form QXj) (see fig. 1) characterized 
by the presence of a marginal fixed point B 0. The corresponding map on the 



torus / is shown in fig. 2, and consists of five complete branches. The domain 
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is accordingly partitioned into five subsets, which we will denote respectively by 
1, 2, 0, 3, 4. Branches 1, 2, 3, 4 have a constant absolute value of the slope A, while in 
region the map takes the usual Manneville and Pomeau formEI x n+ \ = x n + a-x^ 
(7 > 1), leading to intermittency. This partition induces a good symbolic dynamics: 
every possible combination of symbols is physically realized as a trajectory of the 
map. 

To apply the formalism we have to accomplish two things: enumerate all possible 
cycles once the marginal fixed point is pruned away (the effect of the marginal fixed 
point is nonlinear and is probed by the infinity of cycles which accumulate to it@), 
and specify a v and |A P | for each cycle. It is easy to see that the symbolic dynam- 
ics corresponding to prohibiting an infinite repetition of is determined by unre- 
stricted grammar in the following (countable) alphabet {0 fc l, 0*2, n 3, 0*4, 1, 2, 3, 4} 
(k,l,n,i = 1,2,3. . .). As regards winding numbers, if we denote by €j a generic 
letter of the alphabet we have a eie2 ,„ en = J2?=i °<l % where for each letter including 
symbols 1 or 2 a e — 1, while for each letter including symbols 3 or 4 a e — — 1. We 
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now use the piecewise linear approximation of Gaspard and WangU — U to the Man- 



neville Pomeau system to estimate the stability of cycles shadowing the marginal 
fixed point, so that |A eie2 ... e J = lTjLi | A Ej . | , and for the various letters we have 

|A e |=A e = l,2,3,4 

lA^A-A;^ 1 e = k L L = 1,2,3,4 U 

with a = l/( 7 - 1). As each curvatured is zero we get that the zeta function 
is expressed in terms of the fundamental cycles @] (which are determined by the 
letters) and takes the following form 

x 4 cosh (3 ( ~ z k \ 
Co (z, p) = 1 z 1 + 2. 



A V H ka+1 7 

\ fe=l / 

In this way @ reads 

y<t l.-(o+l) 

This leads to normal diffusivity for a > 1, while for other parameter values we get 

o m I t a for a G (0,1) 
a{t) ~\t/\nt for a = l (5) 

in agreement with probabilistic estimates and numerical simulations in 0. We 
remark that our findings are entirely based on metric and topological properties of 
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the system, thus the relation, expressed by fl5|), between anomalous diffusion and 
exponents of power-law stabilities does not rely on any stochastic modelization of 
the system. 

Our theory may also be applied when the diffusion is accelerated EU, that is 
when marginal fixed points appear in the running branches (see fig. 3). We will 
give a detailed treatment of the results elsewhere: here we just quote the results (in 
agreement with flQ[| , except the case a = 1) 



t for a > 2 

t In t for a = 2 

t 3 ~ a for a G (1, 2) (6) 

t 2 1 hit for a = 1 

t 2 for a e (0, 1) 



where a is still 1/(7 — 1). We observe that topological intricacies of the system 
in its laminar phase do not affect the asymptotic estimates (f|) and @, which are 
exclusively determined by the sequence of cycles accumulating to the marginal fixed 
points. It is precisely this robustness with respect to fine details of the dynamics 
which let us claim that our approach may be relevant also in much more complex 
systems, like the Lorentz gas with infinite horizon. This may be appreciated by 
a simple argument: the weight of orbits travelling a time t without collisions in 
this system goes like t _3 E3, corresponding to a = 2 in (|J). From (Bhwe thus get 



<J 2 (t) ~ tint, which indeed seems to reproduce the correct behaviourEH . 

In this paper we have proposed an extension of the periodic orbit theory of 
deterministic diffusion S ^, capable of predicting the asymptotic time behaviour 
of cr 2 (t) by using an appropriate cycle expansion. In particular this technique has 
been tested with classes of one-dimensional maps in which intermittency slows 
or accelerates diffusion: the anomalous diffusion is completely characterized by 
the exponent a ruling power-law stability of cycles shadowing the marginal fixed 
points (a is in turn connected to the intermittency exponent 7). Our results provide 
examples in which cycle expansions can successfully deal with marginal stability: 
this supports the view that this technique may represent a major tool in the analysis 
of generic low-dimensional (both classical and quantum) hamiltonian systems. 

Acknowledgements: We warmly thank Predrag Cvitanovic and Arkady Pikovsky 
for useful comments. 
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Figure captions 

1. Example of a circle map with a marginal fixed point. 

2. Map on the torus associated to the one shown in fig. 1. 

3. / for a model of accelerated diffusion: marginal fixed points are situated in 
the running branches (|<r| = 1). 



6 



